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Abstract 

I present a nonlinear differential equation model (DAM) for the spectrum of Kelvin waves on a 
thin vortex filament. This model preserves the original scaling of the six-wave kinetic equation, its 
direct and inverse cascade solutions, as well as the thermodynamic equilibrium spectra. Further, 
I extend DAM to include the effect of sound radiation by Kelvin waves. I show that, because of 
the phonon radiation, the turbulence spectrum ends at a maximum frequency uj* ~ (e 3 cf /k 16 ) 1 ' 13 
where e is the total energy injection rate, c s is the speed of sound and k is the quantum of 

■ circulation. 
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1 Kelvulence: cascades and spectra. 

i/-) ' Kelvin waves propagating on a thin vortex filament were proposed by Svistunov to be a vehicle for 

■ the turbulent cascades in superfluids near zero temperature [1]. Presently it is a widely accepted 
^ ■ view well supported by the theory and numerical simulations, see e.g. [2-6]. I will refer to the state 

", characterised by random nonlinearly interacting Kelvin waves as "kelvulence" (i.e. Kelvin turbulence). 
^ \ Recently, Kozik and Svistunov [4] used the weak turbulence approach to kelvulence and derived a six- 
\ wave kinetic equation (KE) for the spectrum of weakly nonlinear Kelvin waves. Based on KE, they 
derived a spectrum of waveaction that corresponds to the constant Kolmogorov-like cascade of energy 
?• ■ from small to large wavenumbers, 

•M\ n k ~k~^/\ (1) 

l-H ' Because the number of waves in the leading resonant process is even (i.e. 6), KE conserves not only 
the total energy but also the total waveaction of the system. The systems with two positive conserved 
quantities are known in turbulence to possess a dual cascade behavior. For the Kelvin waves, besides 
the direct energy cascade there also exists an inverse cascade of waveaction, the spectrum for which 
was recently found by Lebedev [7], 

n k ~k- 3 . (2) 



1 



Interestingly, such —3 spectrum was suggested before by Vinen based on a dimensional argument not 
involving the energy flux [2] . Similar argument in water wave turbulence gives famous Phillips spectrum 
which is associated with sharp water crests due to wavebreaking which occurs at large excitation levels. 
By analogy, we could expect that Vinen —3 spectrum should be observed in kelvulence at high excitation 
levels leading to sharp angles due to reconnections. This view is supported by recently reported numerics 
by Vinen, Tsubota and Mitani [5], where it was argued that the —3 exponent arises when the vortex line 
bending angle becomes large (of order one). Kozik and Svistunov [6] later reported a result obtained with 
a refined numerical method which gave a spectrum closer to the —17/5 shape, which is a Kolmogorov-like 
direct cascade of energy dominating the wavebreaking effects at lower turbulence levels. The inverse 
cascade spectrum, although it has the same —3 exponent as the Vinen spectrum, is fundamentally 
different: it corresponds to weak rather than strong turbulence and is more relevant if the main source 
of the Kelvin waves is at small scales. For example, the smallest scales on the vortex filament could be 
generated by reconnections via sharp angles produced by these processes. 



2 Differential equation model 

Differential equation models proved to be very useful for analysis in both weak turbulence [8-10] and 
strong turbulence [12-14]. These equations are constructed in such a way that they preserve the main 
scalings of the original closure (KE in the case of weak turbulence), in particular, its nonlinearity 
degree with respect to the spectrum and its cascade and thermodynamic solutions. For the Kelvin wave 
spectrum these requirements yield, 

where k is the vortex line circulation, C is a dimensionless constant and uo = uo{k) = ^k 2 is the Kelvin 
wave frequency (here, we ignore logarithmic factors). This equation preserves the energy 

E = J u 1/2 ndcu (4) 

and the waveaction 

N= I u-V 2 ndu. (5) 



It is an easy calculation to check that equation (3) has both the direct cascade solution (1) and the 
inverse cascade solution (2). It also has the same thermodynamic Rayleigh- Jeans solutions as the 
original KE, 

T 

n = — — . (6) 

UJ + fl 

where T and fi are constants having a meaning of temperature and the chemical potential respectively. 
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3 Kelvulence radiating sound. 



In contrast with the classical Navier-Stokes flow, there is no viscosity that could dissipate the superfluid 
turbulent cascade at small scales. At low temperatures, friction with the normal component is also 
inefficient and the only dissipative process which can absorb the cascade is radiation of sound by 
moving superfluid vortex filaments [2]. Let introduce the sound dissipation effect into DAM by using 
the classical results of Lighthill about sound produced by classical turbulence [15]. Namely, we will use 
the result that the rate at which the sound energy is generated is proportional to the fourth power of 
the Mach number M = v/c s ~ n 1 ^ 2 /c s , where v is characteristic velocity in turbulence and c s is the 
speed of sound. Further, we will ignore the logarithmic corrections due to the finite vortex core size a. 
In the other words, a should not enter the expression explicitly but only implicitly via c s ~ re/a (which 
in turn enters only via M). Then the rest of the sound dissipation term can be completed uniquely via 
the dimensional argument and the result is 

iu 9/2 n 2 

(fl) radiation i /9 A ' v) 

Let us now examine what effect of the sound radiation on the direct energy cascade. For these purposes 
we can consider even simpler first order DAM which still describes the direct cascade but ignores the 
inverse cascade and the thermodynamic solutions. Such DAM, including the sound radiation term, 
reads: 



Ci wo d , 5 17m C 2 uj 9/2 



n 



2 



>< I2 tt:('^ ;v2 ) - Jr nzr- (») 



re 10 du y ' reV2 c 4 



where C 1 and C 2 are dimensionless constants. The general stationary solution of equation (8) is 



n = A(l-Bu 13 ' 5 ) 1/3 u- 17 ' 10 , (9) 



where 



A = (2/C 1 ) 1 / 5 re 21 / 10 e 1 / 5 , (10) 

B = *-^C 2 k^c-^'\ (11) 

and e is the total energy dissipated in the system per unit time per unit length. 

In the absence of sound radiation, B = 0, we recover the direct cascade spectrum (1). For B > 0, 
the spectrum has the direct cascade shape (1) at low frequencies, u> B~ 5 ' 13 , and it falls to zero at a 
finite frequency 

„ 5- 5 /13 „ (e 3 c 20 /K l 6) l/13_ (12) 

On the other hand, Kelvin waves can only have wavelengths greater than the vortex radius a which in 
terms of the frequency means u> < uj c = c 2 / 're. Thus, expression (9), particularly the finite cut-off at u*, 
will only hold if uj* < uj c or 

e < kc 2 s . (13) 
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This condition can be formulated in terms of the characteristic bending angle of the vortex line a which 
in the direct cascade state is related to e as 



a ~ 



(d 2 /K 3 ) 1/W 



(14) 



where / is characteristic length of the Kelvin waves. Thus, condition uo* < uo c becomes 



a < (I /a) 



1/5 



(15) 



which always holds for a < 1 because / > a. Therefore we conclude that the kelvulence cascade will 
always decay to zero before reaching the maximal allowed frequency of the propagating Kelvin waves. 



In this Letter I presented a differential equation model (DAM) for the Kelvin wave turbulence (kelvu- 
lence) consistent with the scalings of, and having the same set of cascade and thermodynamic solutions 
as, the original six-wave kinetic equation. DAM also exists in reduced versions if some of the flux 
and thermal solutions are not important for a particular problem and can be ignored (in this Letter I 
presented three versions given by the fourth, the second and the first order equations respectively). 

Based on the Lighthill's theory of sound produced by turbulence, I extended DAM to include the 
effect of kelvulence dissipation by sound radiation. I obtained a steady state solution of this model 
corresponding to the direct cascade of wave energy gradually dissipated via radiation and decaying to 
zero at a finite wavenumber u* given by the expression (12). This solution can also be used to predict 
the spectrum of the radiated sound, because the sound energy is just the same as the energy lost by 
kelvulence via radiation. One can also extend this study to scattering of sound by turbulence in systems 
where such sound is generated by externally. 

Finally, in appendix I presented a "warm cascade" solution for kelvulence which is relevant for 
numerical simulations with a finite cut-off frequency. 

I thank Carlo Barenghi, Vladimir Lebedev, Boris Svistunov, Joe Vinen and Grisha Volovik for their 
helpful comments about my model. 

5 Appendix: Warm cascade solutions. 

Above, we ignored the "thermal" component in turbulence which is totally justified because this com- 
ponent does not show up when turbulence is dissipated gradually in the Fourier space, as it is done via 
the sound radiation in our case. On the other hand, it is known in the turbulence theory that sharp 
dissipation or presence of a cut-off frequency lead to reflection of a large portion of the energy flux from 
the smallest scale and, consequently, an pile-up of the spectrum near the smallest scale. This "bottle- 
neck" effect can be described in terms of "warm cascade" solutions in which both an energy flux and a 
thermal components are present. As I showed above, the natural cut-off frequency of Kelvin waves uj c is 



4 Conclusions. 
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not going to be reached by the turbulent cascade because it will always terminate at a frequency < u c 
due to the sound radiation. However, in numerical simulations of kelvulence the cut-off frequency may 
be less than cj* due to the limited numerical resolution. Thus, the warm cascade solutions could be 
relevant for understanding spectra obtained numerically. Such warm cascade solutions were originally 
obtained for classical Navier-Stokes turbulence in [13] and here we will follow a similar approach to find 
such states for kelvulence. To describe both the energy cascade and the thermodynamic component 
(but still ignore the inverse cascade solution) we need to use the second order version of DAM, namely 



where T is an arbitrary constant having the meaning of "temperature". This is the "warm cascade" 
solution which has the pure cascade solution and the thermodynamic energy equipartition as two of its 
limits. Note that the thermodynamic part of this solution describes the spectrum pile-up due to the 
bottleneck phenomenon in numerical simulations due to presence of a cutoff frequency. The relative 
strength of the cascade and the thermodynamic components will be determined by the ratio of the 
incident and the reflected energy fluxes which, in turn, will depend on the particular form of dissipation 
at small scales chosen in numerical simulations. The goal of the true numerical simulation is to chose 
the dissipation function in such a way that the bottleneck effect is minimal or absent. 
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